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Using the non-relativisitc reduction of Coulomb gauge QCD we compute spectrum of the low 
mass hybrid mesons containing a heavy quark-antiquark pair. The gluon degrees of freedom are 
treated in the mean field approximation calibrated to the gluelump spectrum. We discuss the role 
of the non-abelian nature of the QCD Coulomb interaction in the ordering of the spin-parity levels. 
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I. INTRODUCTION 

Gluons are responsible for almost the entire mass of 
light hadrons and contribute significantly to hadron spin. 
Gluons have been probed in hard processes while spec- 
troscopy of gluonic excitations has so far been very lim- 
ited to a few potential hybrid and glueball candidates. In 
the near future, however experiments at JLab, PANDA 
and BESIII are expected to significantly advance our 
knowledge of gluonic excitations through measurements 
of decay channels that couple to exotic quantum numbers 
carried by hybrids. PANDA and BESIII will in particular 
be sensitive to charmonium hybrids. 

Modern studies of gluon dominated hadron spectra 
have been possible through lattice simulations [1-3] while 
in early days of hadron phenomenology various effec- 
tive models have been developed [4-11]. Among these 
the constituent, or quasi-particle approach have had may 
successes in particular in the context of heavy quarkonia. 
For gluonic excitations in presence of heavy quarks the 
Born-Oppenheimer approximation can be employed to 
replace the (fast) gluon field by an effective potential be- 
tween the non-relativistic quarks [12]. These effective po- 
tentials originate from excited gluon configurations and 
can be reliably computed using lattice gauge simulations 
and afterwards used to construct hybrid heavy quarkonia 
[13]. Keeping the quark-antiquark sources static and tak- 
ing their separation to zero results in a spectrum of the so 
called gluelumps which describe the gluon field bound to 
a static color octet source [2] . It is thus natural to expect 
that there should be a close relation between spectrum 
of gluelumps and heavy quarkonium hybrids. This is be- 
cause in the latter the average spacial separation between 
heavy quarks is small, of the order of (r) ~ 0{l/{asm)) 
where, m is the heavy quark mass and as is the running 
couphng evaluated at distance scale 0{{r)). 

Recent lattice simulations of the charmonium spec- 
trum in quenched QCD beyond the Born-Oppenheimer 
approximation also lead to more states then predicted 



by the quarkonium potential models potentially indicat- 
ing presence of gluonic excitations [14]. In the quenched 
approximation and with non-relativistic quarks the con- 
nection between heavy quarkonium spectrum and the po- 
tential model can be made precise. In the context of the 
Wilson loop formalism this has been done in [15-17] 
where potential matrix elements between heavy quarko- 
nium states were related to correlation functions con- 
taining QCD gluon field operators. In [18] the Foldy- 
Wouthousen formalism was used instead and the effec- 
tive non-relativistic QCD Hamiltonian in the Coulomb 
gauge was derived. In this framework all operators con- 
serve the heavy quark quark number and depend on the 
gluon field. The effective QQ potential matrix elements 
in general depend on the distribution of the gluon field in 
the system. In general, gluon distribution in the vacuum 
and in a state containing a QQ pair are different. In [19] 
it was found that at large separations the Coulomb en- 
ergy exceeds that of true static quarkonium state. This 
makes possible to interpret the former as a variational ap- 
proximation to the true QQ QCD state. This difference 
has simple interpretation in the Coulomb gauge [20, 21]. 
The Coulomb energy corresponds to the energy of an ap- 
proximate, variational state with the QQ pair added to 
the QCD vacuum. Gluon-number changing interactions 
in the Coulomb potential couple the QQ source to an 
arbitrary number of transverse gluons which eventually 
produces the string-hke QCD QQ eigenstate whose en- 
ergy is that of the Wilson loop. Extrapolating to small 
QQ separations, however, a reasonable initial approxima- 
tion would be to assume that the gluon field is not signifi- 
cantly disturbed by the heavy quark sources and to build 
the spectrum of gluon excitation using an effective gluon 
Fock space that is independent from the sources. Such a 
spectrum was constructed in [22-28] using a variational 
mean field approximation and describes effective massive 
gluons propagating in a nonperturbative vacuum. We 
will use this Fock space here and combine it with non- 
relativistic quark an antiquark to compute spectrum of 
heavy hybrids. 



In [29] we followed the same approach to compute 
spectrum of gluelumps. By using a single framework of- 
fered by the non-relativistic QCD, we will be able to ex- 
plore the connection between gluelump and hybrid spec- 
tra. In this paper we consider the leading contribution 
to the heavy quarkonium hybrid energies i.e. excluding 
spin dependent-effects. In the following Section wc dis- 
cuss derivation and treatment of the Foldy-Wouthouscn 
Coulomb gauge (FWCG) Hamiltonian and the mean field 
approximation to the QCD vacuum. In Section III 
we summarize calculation of spectrum of ordinary (non- 
hybrid) quarkonia which is used later to benchmark the 
hybrid spectrum which is studied in Section IV. Sum- 
mary and outlook are given in Section V. 



II. FOLDY-WOUTHOUSEN COULOMB 
GAUGE HAMILTONIAN 

The non-relativistic QCD Hamiltonian, HpwCG de- 
scribing interaction between heavy quarks and gluons can 
be constructed from the full QCD Hamiltonian in the 
Coulomb gauge by employing Foldy-Wouthuyscn trans- 
formation [18] 



H -^ HpwcG — THT 



(1) 



Here T = ■ ■ ■ e^S2giSigtSo qj^^ ^j^^ hermitian generators Si 
are of the order {{p) /mY respectively , with m standing 
for mass of the heavy quarks and {p) being the aver- 
age three-momentum of the heavy quarks. The FWCG 
Hamiltonian is thus an effective Hamiltonian valid for 
momenta smaller than the heavy quark mass. The op- 
erators Si eliminate couplings between upper and lower 
components of Dirac spinors to any given order in the 
\/m expansion and lead to a Hamiltonian that conserves 
number of (heavy) quarks. It has the form of a series 
expansion in inverse powers of the heavy quark mass, 



HpwcG — H^i + Hq + Hi + H2 



(2) 



Each term. Hi being of order m* depends on the trans- 
verse gluon and the non-relativistic quark and antiquark 
fields. The transverse gluon degrees of freedom are given 
by the vector potential A°(x) and the conjugated mo- 
mentum n°(x) (V • A" = V • n'' = 0) satisfying 



[A"(x),n^(y)]=<5T(x-y)5 



ab 



where 5T(x-y) = [I- VV/V^](5^(x - y)). The canon- 
ical momentum 11 is the negative of the transverse com- 
ponent of the chromo-electric field. The quark and an- 
tiquark degrees of freedom can be written directly in 
terms of quark creation and annihilation Q\ Q) and anti- 
quark creation and annihilation, {Q, Q^ operators which 
in terms of the original Dirac fields are given by. 



Qa(x) = 



1 + /3 



^(x) 



[Q(x)(-*a2)]A 



^t(x) 



1-/3' 



(4) 



A-l-2 



Here A = 1,2 denotes the z-component of quark spin 
and the Pauli matrix a2 is introduced so that Qx and 
Qx belong the the same SU{2) representation. The mass 
term iJ_i is given by (summation over spin indices is 
implicit) 



H 



dx[gt(x)Q(x)+Qt(x)Q(x) 



(5) 



The Hq term represents the Yang-Mills Hamiltonian cou- 
pled to static quark and antiquark sources. 



Ho = i /dx[^-in''(x)^n''(x)+B''(x).B'^(x)] 
+ ]^ j d^dyJ-^p-{^)K{^,a-y,h)Jp\y). (6) 



Here J — Dot [— V • P] is the determinant of the 
Faddeev-Popov operator; V — Dab = ^afcV + gfacb-^'^ is 
the covariant derivative in the adjoint representation, and 
B is the chromo- magnetic field, B"(x) = V x A"(x) + 
(g/2)/a6cA''(x) X A'=(x). The last term in Eq. (6) rep- 
resent the non-abclian Coulomb gauge interaction be- 
tween color charge densities, p''(x) = Q^{x.)T°'Q{x) — 
0^(x)r*'^Q(x) + /atcA'^(x)n'^(x), and the Coulomb ker- 
nel is given by 



A'(x, a; y, b) = 



^ .(-V^^ ^ 



V-X>' 'V-V 



(x,a:y,6) 



(7) 



(3) The 0(l/m) terms are given by 



^1 = ^ /" dx [Qt(x)D2(x)Q(x) + {Q^Q,T^ T*)] - 1- f d^ [Q^(x).9^ • B''(x)r'^Q(x) - {Q ^ Q,T ^ T*)] 



r 



with D being the covariant derivative in the fundamen- 
tal representation, D = D,y(x) = —i'VSij — gT°jA"-{x.). 



The 0(l/m) Hamiltonian arises from the Z-graph (see 
Fig. 1) and is generated by the first order FW transfor- 
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FIG. 1: Non-relativisitc reduction of heavy quark propagator 
(left side) leads to a kinetic energy term and hyperfine in- 
teraction between quark (chromo) magnetic moment and the 
(chromo)magnetic field. 



mation with Si = {-i/2m) J dxip^x)Pa ■ DV'(x). The 
FW transformation effectively eliminates the quark prop- 
agator in Fig. 1 and shrinks the relative distance between 
the gluon interaction points to zero. This requires that 



J 



the momentum of the quark propagator is smaller then 
m thus, as already mentioned, it is necessary that that 
quark momenta are smaller then the quark mass ( this 
is expected to be the case for low mass heavy quarko- 
nium bound states ). The first term in Hi contains quark 
and anti-quark kinetic energies and 0{l/m) local, spin- 
independent interaction between (anti)quark and gluon. 
The second terms describes 0{l/m) spin-dependent local 
interaction between (anti)quark and gluon. In the lead- 
ing order we drop the 0(l/m) and higher interactions 
terms and keep only kinetic energies of the quarks. In 
Section IV we comment on the results of next to lead- 
ing order calculation, were we add the remaining 0{l/m) 
terms. Finally the second term contains 0(l/m) interac- 
tion that mixes pure QQ states with QQg, hybrid compo- 
nents. These generate 0{l/m^) spin-depdncnt splitting 
which we ignore all together. The final form of Hi which 
is diagonal in the hybrid basis, and is leading order in 
0{l/m) is thus given by 



Hi 



2^ dx [Qt(x)v2Q(x) + gt(x)V2Q(x)] + ^ 



9 
Am 



I dx [gt(x)r''T^A"(x) • A^(x)Q(x) + {Q^Q,T^ T*)] 
J dx [Q^x)T-fabcO- ■ A''(x) X A^(x)Q(x) - {Q ^ Q,T ^ T*)] 



(8) 



with only the first term used in what we define as the 
leading order calculation. 



A. Gluonic degrees of freedom in the mean field 

In the Shrodinger representation the lowest eigenvalue 
of the QCD Hamiltonian ffgcD [A] ^„ [A] = £;„^[A] , 
^o[^] represents transverse gluon field distribution in the 
vacuum. In the Coulomb gauge Fock space construction 
of the excited states (mesons, baryons, hybrids, etc . . . ) 
is motivated by the Gribov-Zwaniziger confinement sce- 
nario. The vacuum state is expected to be dominated by 
field configurations near boundary of the Gribov horizon, 
and gluon excitations which build the hadron spectrum 
correspond to massive quasi-particles of the harmonic ex- 
pansion of the the wave functional near the Gribov hori- 
zon. In practice, in a mean field description [22-28] the 
dominance of the field configurations near the Gribov 
horizon is obtained by choosing the renormalized cou- 
pling .g(A) to be such that the inverse of the Faddeev- 
Popov develops a pole near the zero mode. Consider a 
gaussian ansatz for the vacuum wave functional. 



*o[^] = (A|0)^cxp(^-iy' 



dk 



(27r) 



^A''{k)uj{k)A''{k) 



in the rainbow ladder approximation, show in Fig. 2 the 
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FIG. 2; Diagrammatic representation of the expansion of the 
functional integral for the ghost propagator c.f. Eq. (10) 

expectation value of the inverse Faddecv-Popov operator 



dxe*'^''(*n 



satisfies, {k = k/k), 

1^1 Nc 

d{k) g{A) 2 



V D' 



l*n 



d{k) 



(10) 



'^ (l-k.p)-'^(IP-'^l) 



(2^) 



(9) 



'cj(p)(k-p)2' 

(11) 

The numerical solution to this Dyson equation has been 
extensively studied in [22-28]. The main features of 



its solutions remain unchanged if one uses the angu- 
lar approximation in evaluating the momentum integral 
|p — k| ^ pO{p — k) + k9(k — p) and use the approxima- 
tion uj{k) — > mg{~ constant) for fc ^ and uj{k) — > k for 
k > nig [22]. This approximation enables to obtain an 
analytical solution to Eq. (11) which is given by 



g(A) 



d{k) 



[l + 7l^9^(A)( 



g(A) 



l+(|ff^9"(A)log 



1 1/2 ; '^ ^ ^^g 



T72-, ^ > ''^g 



(12) 



The resulting gap equation involves expectation value 
of the Faddeev-Popov operator whose computation, how- 
ever is not free of from ambiguities. These ambiguities 
arise of UV subtraction which can be absorbed into def- 
inition of Lo [25, 28]. Therefore we choose a more phe- 
nomenological approach. Following our studies of the 
gluelump spectrum [29] we will treat w as a parameter 
of the model to be constrained by the physical spectrum 
rather then by the properties of the vacuum i.e. the gap 
equation. Since the gluon 2-point function is given in 
terms oi lo, by 



with (3h = 8A^c/3 and /3l = 5/3i,/3. It follows that dik) 
can be made renormalization scale invariant by choosing 
51(A) to depend on A in the same way d(k) depends on 
k. This also allows for interpreting dik) as the running 
coupling. For g(A) < gc = 47r(mg//3iA)^/^ the Landau 
pole is avoided and, as g{A) — > gc from below, the 
functional integral of the gluon field in Eq. (10) becomes 
dominated by the configurations near the Gribov hori- 
zon. The parameter uj which fixes the vacuum wave 
functional can be constrained by minimizing the energy 
density (5(*o|ff|*o) =0. 



(*o|A"(k)A''(p)) = (2^)35(k + p) 



Sab 



2uj{k) 



(13) 



it is clear that the wave functional in Eq. (9) leads to a 
gap in single gluon spectrum. The energy of the single 
gluon state defined as 



|l(7) = |k,A,a)=at(k,A,a)|*c 



(14) 



where the single gluon creation (annihilation) operators 
a' ,{a) are defined by (e are hclicity vectors) 



dU 1 

A''(x)= / ., .. , . . [6(k,A)a(k,A,a) + et(-k,A)at(-k,A)]e^'^-^ 



n''(x) 



(27r)3 y2^fcy 
dk juj{k) 



[e(k, A)a(k, A, a) - e^(-k, X)a^-k, A) 



(15) 



is given by [22, 23] 
uj[k) 2 



.^v>,,,|k-kl,i±<!LJ^ 



(27r)3 



2w(fc') 



(16) 

Here Vcl is the negative of the momentum space expec- 
tation value of the Coulomb potential 



VcLik) 



(ixe 



ik-x , 



*o| 



V D 



(V^) 



9 



V D 



1*0) (17) 



which is well approximated by a combination of Coulomb 
and linear potential [22]. The long range part of Vcl 
can potentially lead to and IR divergence in the self en- 
ergy ( integral in Eq. (16)) from the integration region 
k' ~ k. This divergence is however canceled in the 
full bound state equation by similar divergencies aris- 
ing from quark self energies and quark-antiquark inter- 
action. The quadratic UV divergence from k' ^ 00 in the 
self energy integrals is canceled by a gluon mass counter- 
term [22, 23]. In the basis of massive gluons the Foldy- 



Wouthousen Hamiltonian is schematically given by 
Hfwcg ^ 2J Sg(A:)a^(k, A, a)a(k. A, a) -f F(a^a) 

k.A,a 

(18) 
where V represents residual interactions between glu- 
ons and heavy quarks. The gluon mass gap implied by 
Eg(fc) y^ k suggests that a hybrid state in which quarks 
and gluons have residual 3- momenta of the order of Aqcd 
should be well approximated by a state containing a 
pair of non-relativisitc quark and antiquark and a sin- 
gle quasi-gluon. In the infinite heavy quark mass limit, 
separation between the quarks becomes a good quantum 
number and energy spectrum of such static hybrids is de- 
termined by the energy of the single quasi-gluon in pres- 
ence of color octet QQ source. As separation between 
quarks is taken to zero spectrum of static hybrids be- 
comes rotationally invariant and such states are being re- 
ferred to as gluelumps. Lattice simulations show that the 
lowest glue-lump state has spin J, parity P, and charge 
conjugation C, J^'^ = 1"' and the first excited state, 
J^^' = 1 . This unusual parity inversion between vec- 



tor and pscudo- vector states has been explained using the 
Coulomb gauge approach as due to the three-body force 
involving all three particles quark, antiquark and gluon 
originating from the non-abelian Coulomb interaction in 
Eq. (6). A good agreement between our Coulomb gauge 
approach and lattice gluelump and static hybrid spectra 
gives us confidence that the approach may also be ade- 
quate for computation of the spectrum of hybrids with 
dynamical heavy quarks. 



III. ORDINARY QUARKONIA 

Before discussing giuonic excitations we shall fix the 
parameters of the QQ sector by comparing the HpwcG 
Hamiltonian spectrum with that of ordinary charmonia. 



quarkonia. Since low-lying quarkonia are expected to be 
dominated by the QQ component we may represent iV-th 
quarkonium state of spin and z-component J, M parity 
P and charge conjugation C, as 



\JMPCN) 



E 



q^dq 



^':{q)\JMPC;a;q), 



with the QQ state given by 



(19) 



\JMPC;a;q) = ^ I dqxi',!™? (q, 



4™.,n^^-q,™..J0), (20) 



A. Basis and Hamiltonian matrix elements 

The Coulomb gauge Foldy-Wouthuysen Hamiltonian 
should reproduce the known spectrum of ordinary 



and where the spin-orbial wave function x describes cou- 
pling of total quark spin Sq to relative orbital angu- 
lar momentum Lq to total spin of the quarkonium J, 



Lq -|- Sg — J, 



JMPC 
Ami ,^^2 



(q,«) = i[l + C(-l)^'+^']i[l + P(-l)^'+i] Y. {lr,i,;lm2\SqMs){SqMs;LqML\JM)YL,MA^)- (21) 

Ms, Mi 



r 



The Shrodinger equation for the masses of ordinary 
quarkonia is given by 

m 



MMKii) = 



2m+ — 






(22) 

The first term represents quark and antiquark kinetic 
energies and their self-energies, 

dq 



^q ^ -CpVLiO) ^ 'Cf 



(27r) 



rVLiq) (23) 



In position space the expectation value of the Coulomb 
potential given in Eq. (17) leads to a short range part 
VcLiR) -^ Vc{R) oc -l/R for i? ^ and long range 
part, Vcl{R) -^ Vl{R) oc i? for i? — > oo The contribution 
from the short range part is absorbed in the definition of 
the quark mass. The long range part leads to an IR 
divergence which in momentum space and is represented 
by the constant Vl(0) in Eq. (23). This constant, and 
therefore sensitivity to the i? — > oo limit of the position 
space Coulomb interaction is removed by the long-range 
part of the quark-antiaquark potential. 



Vnn 



d2q'PL(q'-q)^CL(|q'-q|) 



(24) 



as it should for color singlet bound states. 

We use the quarkonium spectrum to determine the 
renormalized quark mass. From experimental data we 
find for the the S'-wave {Lq = 0) spin-averaged masses 



M^ = -[Mo- 



3il/i- 



(25) 



M^^ = 3.068 Gov for charmonium and M^. = 9.46 
Gev for bottomonium, respectively. Comparing with re- 
sults of numerical solution to Eq. (22) we extract m^ ~ 
1.16 GeV and to;, = 4.58 GeV, respectively. In Fig. 3 we 
compare the predicted masses of the first radial excita- 
tion of the 5- waves and the ground state (spin-averaged) 
P-waves with experimental data. The observed agree- 
ment is a refiection of the Coulomb gauge kernel correctly 
reproducing the " Coulomb -I- Linear" potential in the re- 
gion of position space covered by the bound quarkonium 
wave functions. 



IV. SPECTRUM OF GLUONIC EXCITATIONS 

As discussed above, gluelumps provide a good bench- 
mark for the spectrum of giuonic excitations. Upon 
neglecting the quark motion (i.e setting Hi = 0) 
(anti)quark positions become good quantum numbers. 



and in particular the spectrum of QQ states with vanish- 
ing relative separation between quark sources is given by 
the energy of gluon cloud. In our approach the low ly- 
ing spectrum is obtained by binding a single quasi gluon 
with kinetic energy given by Eq. (16) to the QQ source 
through the Coulomb kernel given in Eq. (6). More de- 
tails are given in Ref. [20, 21, 29]. The spin of the 
gluelump is then identified with the total angular momen- 
tum Jg of the quasi-gluon, orbiting the static QQ source. 
Once the sources are allowed to move the total spin of 
the resulting hybrid meson is given by the sum of angu- 
lar momenta of all constituents. Since heavy quarks are 
expected to move non-relativisitcaly with average sepa- 
ration of the order of 0{l/asm) it is reasonable to expect 
there will be not much distortion of the gluon wave func- 
tion of the gluelump due to quark motion. Finely, since 
we are ignoring quark-spin dependent interactions it will 
be useful to couple quark spin at the very last to other 
angular momenta before producing a state of good total 
spin J . 



A. Basis and Hamiltonian matrix elements 

From the above discussion it follows that the optimal 
basis of QQg states is such in which the QQ relative angu- 
lar momentum Lq is first coupled to total gluon spin Jg. 
The resulting angular momentum j is then coupled to the 
total quark-antiquark spin Sq to the give the total spin of 
the hybrid state J . In this coupling scheme it is straight- 
forward to take the static quark limit and thus compare 
hybrid and gluelump spectra. Specifically, the A'^-th hy- 
brid state with total spin and its z-axis projection, J, M, 
parity P, charge conjugation C is given by 



(271)3(2^)3 "^'^^ 



(26) 



\JMPCN) = XI / 

a={(J,Sq,Lq,Jq) 

X \JMPC;a;k,q) 



Here the quark-antiquark-gluon, QQg state is given by 
\JMPC;a;k,q)= ^ y'dkdqx^^^?,.(k, q, «) 



mi ,m2 ,(T 






The spin-orbital wave function x describes the (Lq+Jq) + 
Sq coupling discussed above , and cr = ±1 represents 
gluon helicity 



JMPC .]^ ^ ) 
1 



^^^[l + C(-l)^'+^' + l] 



X {JgMg-LqMi\jm){3m-SqMs\JM)Yi^^MM) 



X i?;;:._<,(k)[<5.4+p(-i)-' 



(28) 



Finally, q is the relative momentum between the 
quark-antiquark and k is the momentum of the gluon 
in the hybrid center of mass frame. The spin-orbtal wave 
function is kinematical, i.e. determined by the rotational 
symmetry of HpwcG Hamiltonian, while the radial wave 
function ^^(fc, q) will be determined by diagonalizing the 
Schrodinger equation for the hybrid bound state. Parity 
and charge conjugation are also kinematical and given by 



P 

C 






(29) 



respectively. Here ^ = -|-1 for TM (natural parity) and 
^ = — 1 for TE ( unnatural parity ) gluon state which 
correspond to a |cr = +1) +i^|(t = —1) combination of 
gluon helicity states. We note that, as expected, both 
P and C are a product of QQ and gluelump parity and 
charge conjugation given by 



Pq = (-l)^'+\ Pg^a-^V' 
Cq = (-1)^'+'^% Cg^-1. 



(30) 



(27) 



The Schrodinger equation for the radial wave function 
has the form of 



2m+ — + -— + S„(fc) 
m 4™ 



*a(fc,9) 






(31) 



r 



In leading order the fist term represents kinetic energies 
and self-energies of the gluon, given by Eq. (16), and of 
quarks given by Eq. (23). As in the case of Eq. (22) the 
IR divergencies in the self energies are canceled by the 
potential matrix elements. The first integral in Eq. (31) 
is the QQ, interaction. It is independent of the gluon de- 
grees of freedom and repulsive in for the QQ pair in the 
color octet. Finally, the last term on the left hand side 
of Eq. (31) represents collectively the two body attrac- 



J 



tive interactions between the quark and the gluon and 
the antiquark and the gluon and the irreducible three- 
body interaction linking all three constituents - quark, 
antiquark and gluon. 



^QQ9 



y2_ 
QQg 



y3_ 
QQg- 



(32) 



The two-body interaction is given by 



^ooo(^'9'";'=''9':a) 



\ <jj{k) ujjk') 

uj{k') ^ V Lj{k) 



(33) 



7ni ,1712.(7, a' 

X [S{q + k/2 - q' - k'/2) + 5(-q + k/2 + q' - k72)] VclO^ - k')e-k',a' • elk 

I 



The three body interaction emerges from the non-abelian 
Coulomb kernel in Eq. (7). The general expression for the 
three body interaction, V^^ is quite complicated, and 
we discuss special cases in the Appendix. The matrix 
elements of next to leading order interactions are also 
given in the Appendix. 



B. Hybrid vs gluelump spectra 

The particular choice of the coupling of spin and orbital 
degrees of freedom in Eq. (28) enables to make a direct 
contact with the gluelump spectrum. In that case the 
orbital part of the QQ wave function is in the S'-wavc and 
the radial wave function is proportional to a (5-function in 
coordinate space. We take the radial part as a normalized 
gaussian, 



*g(9) 



7r4 



27r 



(27r)3 



^Ui)\' 



(34) 

The leading order Hamiltonian truncated to the S -wave 
quark-antiquark orbital is given in the Appendix. It is 
straightforward to verify that removing the quark mass 
term and taking the limit ag ^ oo ( after m ^ oo ) in 
Eq. (31), which corresponds to setting the relative quark 
separation to zero, the gluelump Schrodindgcr equation 






from [29] is reproduced. 

With dynamical quarks the quark-gluon interaction 
mixes the quark orbital angular momentum with gluon 
spin while in leading order it preserves quark spin. This 
is represented by a mixing between a and a' indices. To 
leading order, however, quark spin is conserved, Sq 
and then parity and charge conservation imply no mix- 
ing between TE and TM gluons. As long as the quark 
motion does not distort much the gluon distribution of 
a gluelump one would expect mixing between various 
{Lq,Jg) states to be small. Spin splitting is also ex- 
pected to be small, of the order of 0{l/m). In other 
words spectrum of low hybrids are expected to be ap- 
proximately classified by the product of gluelump and 
QQ quantum numbers. In the following we will com- 
pute the hybrid spectrum under this approximation. Do 
not see much point in going beyond this approximation. 
This because we expect uncertainty in the spectrum as- 
sociated with our mean field treatment of gluonic degrees 
of freedom to be of similar order of magnitude. This will 
be shown in Section IV C were we study dependence of 
the hybrid spectrum on the mean field parameter, u; {cf. 
Eq. (9)). We also restrict our study to lowest multiplets, 
since bare bound state spectrum is expected to be renor- 
malized above open decay channels. 

The two lowest {Lq,Jg) multiples , {Lq,Jq) = 
(0, 1), (1, 1) lead to the following set of states. For Jg — 1 



there two possible gluclumps, the ground state with 
J^'-^ = 1"' and gluon in the TE mode and the first ex- 
cited gluelump state with J^'-^ = 1^^ and the TM gluon. 
Coupling the TE gluon with the color octet QQ state 
in L„ = 0, S-wave orbital leads to a hybrid state with 



the intermediate angular momentum j — ^^^ 



L„ + Jg = 1. 



After adding the quark spin Sq = 0, 1, in absence of 
spin-dependent interaction and weak mixing between dif- 
ferent quark orbitals. we obtain four low lying hybrids 
with quantum numbers, J 



PC 



r 



for Sq = Q and 



J^^ = 0-+, 1"+, 2"+ for Sq = 1. It is worth noting that 
hybrid with exotic quantum numbers 1 ^ appears in this 
lowest multiplet and is predicted to have the QQ pair in 
spin-1. The TM gluon coupled to the octet QQ state 
would lead to four hybrids with J^'~^ 
and JP^ 



0++,l++,2++ for Sq = 1. 



= 1+- for Sq^O 
The L„ 



.q-lQQ 

multiplet, is expected to be at a similar energy if orbital 
excitations of the quarks is energetically comparable to 
orbital excitation of the quasi-gluon. 



C. Numerical Results 

Solving the radial Shrodinger equation, (31), even for 
fixed Lq , Jg is quite challenging numerically. We proceed 
with the following numerical approximation. First we 
use a single harmonic oscillator quark wave function to 
solve for the low-lying gluon energy in the background of 
either 5* {Lq = 0) or P {Lq — 1) wave quarks. The size of 
the quark wave function is determined by minimizing the 
total hybrid energy. In the next step and to validate sin- 
gle harmonic oscillator approximation for the quark wave 
function , we use the solution for the gluon wave function 
from the first step and solve for the quark wave function. 
We find variation in the hybrid energies to be less then 
a percent. Taking advantage of the gaussian quark wave 
function, we can introduce an effective potential for the 
gluon field and, as discussed earlier make the connection 
between hybrid and gluelump states. In this approxima- 
tion hybrid can be viewed as a single quasi-gluon mov- 
ing in an effective potential originating from the slow 
moving quarks. In the Appendix we give explicit formu- 
las for the gluon interactions involving S'-wave quarks. 

In Fig. 3 full boxes show our results for the spin- 
averaged masses of ordinary quarkonia. As discussed in 
Sec. Ill the ^ and 1 states were used to fix the quark 
quark mass. We find the first radial 5-wave excitation 
at Mqi-+ =3.82 GeV and the center of mass of the P- 
waves at Mi+- = Mq++ = Mi++ = M2++ = 3.53 GeV, 
which compare favorably to the experimental values as 
seen from Table I. 

As discussed in [29] we determine gluon mean field 
parameter u;{k) which appears in the ansatz for the vac- 
uum, Eq. (9) by fitting the gluelump spectrum. Using 
ui{k) specified by model-3 in [29] for the masses of low- 
est two gluelump states J^'^ = 1^ and 1 , gives 
0.89 GeV and 1.29 GeV compared to lattice results of 
0.87(15) Gev and 1.25(16) Gev, respectively. In compar- 
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FIG. 3: Charmonium (solid boxes) and charmonium hybrid 
spectrum compared with data (where available) or lattice 
computations. Single dashed boxes are the ccg hybrids dom- 
inated by the P-wave quarks, all other have the QQ pair in 
the relative S-wave orbital. 
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FIG. 4: Same as Fig. 3 for bottomonium. 



ison, the original vacuum model of Rcf. [22, 23], referred 
to in [29] as model- 1 results in higher gluelump masses 
of 2.05 GeV and 2.19 GeV, respectively. 

We can use the gluelump and ordinary QQ masses to 
obtain a simple estimate for the hybrid spectrum. For hy- 
brids with the TE gluon i. e. built on top of the ground 
state, 1"' gluelump and the S'-wave QQ state we expect 
the exotic 1 '' charmonium hybrid at M2-+ ~ rrig + 
2mQ + mQ{v^) ~0.87(2.05)-f2.3-h0.5= 3.67(4.85) GeV 
for model-3 (model-1) were we used (w^) ^ 0.4 ~ 0.5. 



TABLE I: QQ charmonium spectrum 



J^^' 


Exp. GeV 


This work GeV] 


J^^' 


Exp GeV 


This work GeV 


0-+ 


2.980(1) 


3.07 


0'-+ 


3.637(4) 


3.82 


1-- 


3.097(0) 


3.07 


1'— 


3.686(0) 


3.82 


1+- 


3.526(0) 


3.53 


l'+- 


- 


4.16 


0++ 


3.415(0) 


3.53 


0'++ 


- 


4.16 


1++ 


3.511(0) 


3.53 


1'++ 


- 


4.16 


2++ 


3.556(0) 


3.53 


2'++ 


3.929(5) 


4.16 




FIG. 5; Vacuum wave functional parameter tj(fc), model-l 
and model-3 were discussed in [29]. The dotted line is used 
here and give the best fit to the lowest gluelump and 1 '" 
hybrid. 



ccg 

S-Wave 





II 


II 






1 


1 1 










to»J J-bodv 




- 










^^H 2-body 




- 




















- 










- 


- 


t^ 


rAAA^W^AV.^W^j 


^ 
""i^ 


E 




- 


III 










^ 




















II 


II 




1 1 1 1 





f^ 0" l" 2 T 0"* r* 2* 3"* 0*" 1*" 2*' 3*" O" l" 2" 3*' 



The recent lattice calculation of [14] palce the 1 ^ 
charmonium exotic at 4.33 GeV, i.e. between our two 
models of the vacuum wave functional. Thus finally we 
vary uiik) to obtain masses closest to both the ground 
state gluelump and the 1 '" hybrid. The correspond- 
ing w is shown in Fig. 5 for which we find the 1 ^ 
hybrid at 4.47 GeV and the two lowest gluelumps at 
Mi+- = l.72Gev and Mi— = 2.l2Gev. 

The final results for the charmonium and bottomonium 
spectrum are shown in Fig. 3 and Fig. 4 respectively. As 
discussed earlier, if mixing between (Lg, Jg) multiplets is 
excluded and with hyperfine interactions turned off we 
find 4 QQg degenerate states for Lg = and 12 degen- 
erate states for Lg = 1. We have computed the effects of 
0{l/m) interactions which include spin dependent terms 
and find them to contribute at the level of 30 — 50 MeV 
to the charmonium states which is certainly below the 
accuracy of the variations in the mean field wave func- 
tional. Our results are compared to experimental data 
and/or lattice calculations where available. The particu- 
lar degeneracy pattern we find seems to suggest that the 
PC = — h triplet of states found in lattice simulations 
at 4.3 GeV does indeed contain the 1 ^ exotic [30]. In 
Tables II, III we give our predictions for the masses of 
the ccg hybrids and the possible identification with the 



FIG. 6: The effect of the irreducible three body potential on 
the charmonium hybrid spectrum with S-wave quarks. Sold 
lines represent the spectrum form Eq. (31) with V''^ in Eq. (32) 
removed. 



lattice states from [14]. 

Finally in Fig. 6 we show the effect of the irreducible 
three body interaction, Vqq , that is germane to the 
Coulomb gauge. It is this interaction that is respon- 
sible for producing the inverted parity ordering of the 
gluelump spectra, i.e producing the l"* gluelump be- 
low the 1 . It leads to analogous parity-inversion 
in the QQg quarkonium spectra, i.e. it leads to the 
1--, 0"+, 1"+, 2-+ multiplet below 1+^,0++, 1++, 2++. 



V. SUMMARY AND OUTLOOK 

We studied spectrum of heavy quarkonium hybrids in 
Coulomb gauge QCD with gluon degrees of freedom in 
the mean field approximation. In general, we have found 
a reasonable agreement with lattice data. In particular 
we find the 1 '" exotic charmonium state at 4.47 GeV. 
Our predictions seem to be systematically higher com- 
pared to lattice which is most likely a reflection of a de- 
ficiency of the gaussian, mean filed, vacuum wave func- 
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TABLE II; Charmonium hybrids with S-wa.ve QQ pair compared with lattice results were avilable. The second column shows 
masses all states in a given multiplet (third column) from the Coulomb gauge calculation. In the forth column we give the 
lattice results that are closest to our predictions. The states marked by [?] (if exist) have not been resolved on the lattice. 
Lattice computations find lower mass states for these quantum numbers, it is possible that these could be interpreted as QQ 
states. The exception is the 3 '' for which no lattice state was found. 



■Jg 


This work [GeV] 


J^^ 


Lattice [14 [GeV] 


1+ 


4.476 


0- + ,l-+,2-+, 1-- 


4.291(48),4.327(36),4.376(24), ? 


1" 


4.762 


l+-,2++,[0++,l++ 


4.521(48),4.508(48), [?,?] 


2+ 


5.144 


l" + ,[2--,2-+,3"+ 


4.696(103), ?,?,? 


2- 


5.065 


2+-, l++,2++,3++ 


4.733(42), ?,?,? 



TABLE III: Same as in Table II for QQ in the P-wave 





This work [Gev] 


J^^ 


Lattice [14] [Gev] 


1" 


4.886 


0" + ,l"+,2-+, 7 more 


4.291(48),4.327(36),4.376(24), ? 


1+ 


4.692 


0+^,1+^,2++, 7 more 


4.521(54),4.521(48),4.508(48), ? 


2" 


5.221 


1" + , [11 more] 


4.696(103), [?] 


2+ 


5.276 


2+", [11 more] 


4.733(42), [?] 



tional. This needs to be improved in particular by the 
effects of vortex field configurations, wildly believed to be 
revenant for confinement. The framework however seem 
to be well suited for further investigations of quarkonium 
and hybrid structure, including mixing and transitions. 
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APPENDIX A: HAMILTONIAN MATRIX 
ELEMENTS 



The matrix elements of Vqq and Vqq in Eq. (32) 

are computed using gaussian wave functionals for the QQ 
radial wave function. 



q'^dq q'^dq' L,, w,2,3 



^' «' ^/\m^if^'\ 



(2.)3 (2^)3 ^Q ('^)^Q0.(^' g' "' fe ' '^ ' " )^Q (^ ) 

(Al) 



In particular for the S'-wave, given by Eq. (34) wc find 



J 



yQQ,ik..k') = Sj,j^N, 



I ujjk) L0{k') 

oj{k') ^ Y uj{k) 



^)T.\n+''i-mv-i-m' 



^ _ik-k'[2 

X (27r) / dk-k'P,(k-k')VcL(|k-k'|)e ''"Q . 



(A2) 



and 



yQQ.i^y) = -x/'^'^^'^'T^t^ ^ X™^,r.^.(k,ak4^1S^'(k',«') 



,m2,cr,CT' 



X (47r) 



j3„ _^p+k+k^ 

"^^ i^(^)(|k' + p|,|k + p|,|p|)[2e ^^^ 



k-k'|- 



(27r 



]p-e_fc',^'P-e_fc,o 



(A3) 



where [29], K'-^\k,p,q) — VcLik)Vc{p)Vc{q) + parity and charge conjugation satisfying P = (—1)' i.e. 
Vc{k)VcL{p)Vc{q) + Vcik)Vcip)VcLiq). With hybrid 



11 



the possible values of I are constrained by parity of a 
state, C = (—1)'^^"'^ where is S* = 0, 1 is the spin of the 



QQ pair. To order 0{l/m) the other two interactions in 
Eq. (8) give 



V' 






2Nc ^JuJyUk 



X / dk-k'Pj(k-k')(47r)2a(|k-k'|)e ''"Q 



for spin-idependent contact interaction and, 



V,%{k,k') = ^_^^[l + P(-l/](-l)i+'+.^(25+l)(2; + l)(ll;Z0|J,l) 



Am ^UJkUJk ; 
5 1 1 



2 2 

5* 1 



1 t t 11 i "§ "^' f / rfk.k'P,(k-k')(4^)^«(|k-k'|)e 



r 



for spin-dependent contact interaction. Here 

AnZ 
a(p) 






(A4) 



with Z = 5.94, c = 40.68, Aqcd = 0.25 GeV arising from 
the fit to the Coulomb part of the QCD potential [29]. 
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